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EXTREMELY METAL-POOR STARS AND A HIERARCHICAL CHEMICAL EVOLUTION MODEL 
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Draft version January 20, 2013 

ABSTRACT 

Early phases of the chemical evolution and formation history of extremely metal poor (EMP) stars 
are investigated using hierarchical galaxy formation models. We build a merger tree of the Galaxy 
according to the extended Press-Schechter theory. We follow the chemical evolution along the tree, 
and compare the model results to the metallicity distribution function (MDF) and abundance ratio 
distribution of the Milky Way halo. We adopt three different initial mass functions (IMFs). In a 
previous studies, we argue that typical mass of EMP stars should be high-mass(^ 10 M©) based 
on studies of binary origin carbon-rich EMP stars. In this study, we show that only the high-mass 
IMF can explain a observed small number of EMP stars. For relative element abundances, the high- 
mass IMF and the Salpeter IMF predict similar distributions. We also investigate dependence on 
nuclcosynthetic yields of supernovae (SNe). The theoretical SN yields by Kobayashi et al.(2006) and 
Chiefh & Limonge (2004) show reasonable agreement with observations for a-elements. Our model 
predicts significant scatter of element abundances at [Fe/H] < —3. Best fit yields for one zone chemical 
evolution model by Francois et al.(2004) well reproduces the trend of the typical abundances of EMP 
stars but our model with their yield predicts much larger scatter of abundances than the observations. 
The model with hypcrnovae predicts Zn abundance in agreement with observations but other models 
predict lower [Zn/Fe]. Ejecta from the hypcrnovae with large explosion energy is mixed in large mass 
and decreases scatter of the element abundances. 

Subject headings: stars: abundances - stars: Population II - Galaxy: formation - Galaxy: evolution 



1. INTRODUCTION 

Extremely metal-poor (EMP, [Fe/H] < -2.5 in this 
paper) stars arc stars formed in the early universe in 
terms of chemical evolution. They are thought to have 
been formed at high rcdshift but are still shining with 
the glow of nuclear burning in the Local Group. Re- 
cent large-scaled surveys have identified hundreds of the 
EMP stars in the Milky Way halo (HK survey. Beers et 
al. 1992: Hamburg/ESO [HES] survey, Christlieb et al. 
2003, 2008: Sloan Extension for Galactic Understanding 
and Exploration [SEGUE], Yanny et al. 2009). Element 
abundances of these stars are revealed by follow-up spec- 
troscopic observations. They provide a means of probing 
the earliest phases of the evolution of the Milky Way and 
supernovae (SNe) in the early universe. 

In this paper, we refer to stars with [Fe/H] < —2.5 
as EMP stars, al though "E MP star" is usually used 
for [Fe/H] < - 3 (iBeers et a l. 2005). As shown' in our 
previous study (jKomiva et al.ll2010L hereafter Paper I), 
stars with [Fe/H] < —2.5 show some observational and 
theoretical peculiarities distinguishing them from more 
metal rich Population II stars (see Section 2.1 of Paper 
I). Esp ecially, in previous studies (jKomiva et al.l [20071 
l2009aD . we show that the initial mass function (IMF) of 
stars with [Fe/H] < —2.5 should differ from metal rich 
stars. Typical mass of the EMP stars is ^ 10 M© and 
present EMP stars in the Milky Way halo are the low 
mass minorities. We refer to the mother stellar popu- 
lation with [Fe/H] < —2.5 as EMP population and low 
mass survivors with nuclear burning now as EMP sur- 
vivors. Stars with [Fc/H] < —5 arc referred to as hy- 
per metal-poor (HMP) stars. 2 HMP stars detected in 
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the Milky Way halo are the most metal deficient objects 
observed yet (HE1327-2326, Frebel et al. 2005: HE0107- 
5240, Christlieb et al. 2002). 

Formation environment of the EMP population stars 
are thought to differ from present galaxies. EMP stars 
are formed in the process of galaxy formation in the early 
universe. In the A cold dark matter (CDM) universe, 
large galaxies like the Milky Way are formed through 
merger of smaller galaxies as building blocks. According 
to the hierarchical structure formation scenario, a stel- 
lar halo is an aggregation of stars formed in t he many 
small galaxies ()Searle fc Zinnlll978l:lHelmill2008[ ). Earlier 
theoretical studies show that the first st ars are formed 
in ve r y low mass halos w ith ~ 10^ M© (jTegmark et al.l 
[19971 iNishi fc Susal[T999h and host galaxi es of the sec- 
ond generation of stars are also small (jRicotti et al.l 
[200l lWise fc Abe]|[2Ci08[ ). Chemical abundances of these 
building blocks can differ from each other. Unlike metal 
rich stars, metals in EMP stars are synthesized by only 
one or a few precursory SN(e). Element abundances of 
the EMP stars can reflect individual characteristics of 
the precursory SN(e) and their host galaxies. A semi- 
analytic hierarchical approach can provide a framework 
within which to study the earliest phases of the chemical 
evolution and the formation history of the EMP stars. 

One important point at issue for the earliest phases 
of chemical evolution i s a possible difference in the 
IMF of EMP stars (e.g. lAbia et al.l [20011: iKomiva et al] 
I2007D . Theoretically, typical mass of stars is large in ex- 
tremely metal-poor environment and / or in the early uni- 
verse. Numerical simulations show that Population HI 
stars without metal are very massive (e.g. iBromm et al.l 
119991: lAbel et al.l |2002[) . For EMP stars with a little 
metal, the IMF can be different from nearby stars, too 
(jOmukai et al.l [20051 ). Existence of the low mass EMP 
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survivors in the Galactic halo proves that some low mass 
EMP stars can be formed, but typical mass of the EMP 
popula tion stars can b e more massi ve than Population I 
stars. iKomiva et all (|2007l . l2009a[) give constraints on 
the IMF of EMP stars from statistics of observed EMP 
survivors. It is known that large fraction (~ 20%) of 
the EMP survivors comprise carbon enhanced stars re- 
ferred to as Carbon Enhanced Metal-Poor (CEMP) stars. 
More than half of them show large s-process element en- 
hancement. Abundance anomalies of the s-process ele- 
ment enhanced CEMP stars (CEMP-s stars) are due to 
binary mass transfer. Intermediate massive EMP stars 
with 0.8 — 3 Mq synthesize carbon and s-process elements 
in the asymptotic giant branch (AGB) phase, and pollute 
their companions to make them CEMP-s stars. We pro- 
pose that CEMP stars without s-process enhancement 
(CEMP-nos stars) are also formed through binary mass 
transfer but from more massive primaries with 4 — 6 Mq . 
From observational statistics of the CEMP-s and CEMP- 
nos stars, we give constraints on mass distribution of 
primary stars of EMP binaries an d conclude that typ - 
ical mass of EMP st ars is large (jKomiva et al.l 120071 ). 
iKomiva et al.l ()2009al ) discuss the constraints of the IMF 
from CEMP stars again in detail and give an additional 
constraint from the total number of EMP survivors. The 
number of EMP survivors in the Galactic halo is very 
small and it indicates that the fraction of low mass sur- 
vivors among the EMP population is small. As a result, 
a lognormal IMF with medium mass, A/md ~ 10 M0, 
and dispersion, Am 0.4, can satisfy all the constraints. 
iLucatello et al.l ()2005l) also give a constraint on the IMF 
of EMP stars from statistics of CEMP-s stars and ar- 
gue that the typical mass of EMP stars is slightly higher 
than more metal rich stars. However, because they do 
not take account of CEMP-nos stars and assume a stellar 
evolution model different from ours, they conclude lower 
typical mass, Mmd = 0.79 M0. 

In paper I, we build a merger tree of the Galaxy and 
compute the enrichment history of iron abundance along 
the tree. The hig h mass IMF with M^d = 10 A/© 
(jKomiva et al.|[2007 ) is adopted in the computations. We 
also discussed origin of HMP stars considering the effect 
of surface pollution by accretion of metal enriched inter- 
stellar matter. In this paper, we investigate chemical 
evolution of several elements with detailed theoretical 
nucleosynthetic yields of metal deficient massive stars. 
Model results with high and low mass IMFs are com- 
pared with compiled observational data. To deal with 
individual characteristics of the EMP stars, all the indi- 
vidual EMP population stars are registered in our com- 
putations. We discuss not only averaged abundances but 
al so dispersion of th em . 

iTumlinsonI (|2006f ) presents a semianalytic hierarc hical 
model for the Galactic halo. ISalvadori et al.l (|2006l ) also 
calculate a hierarchical model with gas blowout from ha- 
los by bursty star formation taken into account, but they 
do not deal with individual stars and do not investi- 
gate diversity of the element abundances. These previous 
studies assume the Salpeter IMF for EMP stars. Addi- 
tionally, these previous studies with hierarchical models 
investigate only iron abundance. We calculate forma- 
tion and evolution of stellar halo with different IMFs 
and compare the predicted metallicity distribution func- 
tions (MDFs) and abundance ratio distributions. We use 



four different sets of core-collapse SN yi elds calculated 
for me tal-poor stars. lArgast et al.l (|2000[ ) and iKarlssoiJ 
(|2005l ) investigate element abundances of EMP stars us- 
ing inhomogeneous chemical evolution models. However, 
they do not take account of the merging history of the 
Galaxy and stars are assumed to be formed randomly in 
space. They also do not investigate the IMF dependence. 

This paper is organized as follows. Computation 
method and assumptions appear in the next section. Es- 
pecially, assumptions about IMFs and SN yields are de- 
scribed in Sections l2.2l and l273l respectively. In Section|3l 
observational sample for comparison is described. We 
give results in Section |4] and conclude the paper in Sec- 
tion [5] 

2. COMPUTATION METHOD 

2.1. A Hierarchical Chemical Evolution Model for EMP 

stars 

A hierarchical chemical evolution model for EMP stars 
is developed in Paper I. We built a merger tree semi- 
analyt ically using the method of ISomerville fc KolattI 
( 19991) based on the extended Press-Schechter formalism 
( Lacev fc Cold [r993| ) . Along the merger tree, chemical 
enrichment and formation history of EMP stars are cal- 
culated. In this paper, abundances of O, Na, Mg, Si, Cr, 
Fe, and Zn are computed and compared to observations 
of metal poor stars. 

Stars are assumed to be formed in halos with virial 
temperature, Tvir, higher than lO^K. Star formation ef- 
ficiency assumed to be constant and determined to give 
[Fe/H] = at z = (2.1 x 10"" - 1.1 x lO'^Vyr for 
the following computations). To investigate diversity of 
element abundances, all the individual EMP population 
stars are registered in our computations. Mass of each 
EMP star is specified randomly according to the statis- 
tical weight with the IMF. Half of the all stellar systems 
are taken to be binaries and a flat mass ratio distribution 
is assumed. Adopted IMFs are described in Section [2^ in 
detail. We set nucleosynthesis yields of each massive star 
as a function of it's initial mass and metallicity. Assump- 
tions about yields are described in Section 12.31 Each 
minihalo is assumed to be chemically homogeneous. 

We use the same assumptions as Model P in Paper I 
about radiative and dynamical feedback from massive 
stars. Massive stars ionize the ambient matter. At 
z < 20, Lyman- Werner background radiation inhibits 
star formation in mini-halos those are not pre-ionized. 
Energetic SN explosions blowout gas in their host mini- 
halos when their effective kinetic energy, eE\^, is larger 
than binding energy, E^bin, of the gas of their host halos. 
Gas and metal ejected from minihalos are mixed imme- 
diately and homogeneously throughout the intergalactic 
medium (IGM). Assumed explosion energies, Ek, of SNe 
are described in Sec. 12.31 From larger halos, a little frac- 
tion, r]eEk/ Etiin, of metal ejected by SNe is assumed to 
go to IGM, where e is the fraction of supernovae explo- 
sion energy converted into kinetic energy, and 77 is the 
fraction of the input kinetic energy that is retained by 
gas that escapes their host halo by wind. We assume 
e = 0.1 and 77 = 0.1, but results in this paper arc almost 
independent from these parameters. Because of reioniza- 
tion, mini-halos with Tvir < 10'*K cannot accumulate gas 
at z < 10. 

Computations using different IMFs and different nucle- 
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Fig. 1. — Initial mass functions (IMFs) adopted in our compu- 
tations. The solid red , dashed green and dott ed b lue lines denote 
IMFs of lKomiva et al.l ||2007|' ) : [Chabrieil ||2003|) and lLucatello et al.l 
||2005I '|. respectively. 

osynthetic yields are presented. Adopted IMFs and SN 
yields are described in following subsections and summa- 
rized in Table 1. 

2.2. IMF 

As stated in Sec.l, in our earlier studies, we give con- 
straints on the IMF of the EMP population stars from 
statistics of EMP survivors. We assume lognormal form. 
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and conclude that the high mass IMF with M„^d ~ 
1OM0,Am = 0.4 is in agreement with all the observa- 
tional features. We use this high mass IMF as the fiducial 
one. 

Based on recent observations. IChabrieil ()2003f ) presents 
an IMF of the Galactic spheroid; a lognormal IMF with 
Mmd = 0.22 Mq, Am = 0.33 for stars with m < 0.7 Mq, 
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for higher mass stars. For compar- 



ison, computation with this standard low mass IMF is 
also presented (Model CK). Since this IMF is the same 
as the Salpeter IMF for high mass and intermediate mas- 
sive stars, chemical evolution under this IMF is almost 
the same as the Salpeter IMF. 

The IMF derived by iLucatello erall (I2005D is also 
tested (Model LK). They argue a lognormal IMF peaked 
with slightly higher mass, Mmd = 0.79 M©, Am = 0.51, 
based on the statistics of the CEMP-s stars. 

For simplicity, we use the same IMF for all the stellar 
populations in each computation. But we note that the 
high mass IMF is derived from statistics of EMP sur- 
vivors. In the Galaxy today, more metal rich stars are 
formed under a low mass IMF. 

The IMF of the first stars can be different from EMP 
stars. In Paper I, we have discussed typical mass of 
the first stars from viewpoint of the number of surviv- 
ing Population III stars. Observational scarcity of stars 
with [Fe/H] < —4 suggests that the typical mass of the 
first stars in each mini-halo is higher than the subse- 
quent generations of stars. In this paper, A/md = 40 Mq 
is assumed for local first stars without SN progenitor in 
their host mini-halos. The number of the predicted HMP 
stars becomes comparable with observations under this 
assumption. 

Many numerical studies for the Population III star 
formation argue that one very massive single star 



(or binary, iTurk et al.l (|2009[ )) with m > 100 Mq 
is forme d at the center of the primordial min i- 
halo fc.g. lO'Shea fc Normanl[2007t lYoshida^raII [2008^. 
On the other hand, some studies show that zero- 
metallicity gas can fragment into multiple pieces 

and some low-mass stars a re als o formed withou t 

metal (|Nakamura fc Umemural 120011: IClark et al.l |2008i ). 
Also in the latter studies, however, the typical mass 
of Populatioii III s tars is higher than nearby stars. 
lYoshida et al.l ()2007() shows that stars with ~ 40 Mq are 
formed in photo-ionized halo without metal. When a 
Population III star with ~ 200 Mq is formed, it explodes 
as pair instability SN(PISN). We discuss the effect of the 
very massive first stars and PISNe in Appendix. 

2.3. Supernova Yields 

In this study, stars with 10 — 50 Mq are assumed to 
explode as Type II SNe (SNe II). In Paper I, we simply 
assume that any SNe II eject 0.07 A/q of iron, but in 
this study, we consider dependence of yields on stellar 
initial mass and initial metallicity. For SNe II, we adopt 
four sets of yields. Figure [2] summarizes yields of SNe II 
with Z = against initial mass, by theoretical models 
adopted in this study. 

Theoretical yields by iKobavashi et al.l (|2006f ) are used 
as the fiducial one (Model KK). They give yields of mas- 
sive stars with various initial mass and initial metallicity 
by numerical calculations of stellar evolution and explo- 
sive nucleosynthesis at SNe II. They present yields of 
hyp ernovae with larger ex plosion energy, too. Follow- 
ing |Kob^aihreF^ (|2006f ). in this paper, 50% of stars 
with TO > 20 Mq are assumed to explode as hypernovae. 
Explosion energy of the SNe II are assumed to be Ey^ = 
lO^ierg for normal SNe and i?k = 10^^ x (m/ Mq - 10)erg 
for hypernovae. The most prominent nucleosynthetic fea- 
ture of the hyperno ya is a large Zn y i eld. W e note that, 
in computation by IKobavashi et al.l ()2006( ). iron yields 
for normal SNe are tuned to 0.07 A^© by choosing the 
"mass cut" parameter. In such computations for SN nu- 
cleosynthesis, the location of the boundary between the 
part of the star that eventually collapses to a compact 
object and that which is expelled outward is a free pa- 
rameter and referred to as mass cut. An amount of the 
ejected iron strongly depends on this parameter. For hy- 
pernovae, parameters involved in the mixing and fallback 
are determined to give [0/Fe] = 0.5. We also present 
a chemical evolution model without hypernova (Model 
KKn) to show the contributi on from hypernovae . 

The theoretical yields by IWooslev fc Weaver! ()1995D 
are also adopted (Model KW). We use the explosion 
m odels labeled as 12A-22A 25B, 30C, 35C and 40C 
in IWooslev fc Weaverl ()1995[ ). For explosion energy, re- 
sultant kinematic energy of their computations is used 
( 1.1 - 3.01 X lO^ ^er g). 

iFrancois et al.l (l2004f) modify the yields by 
IWooslev fc Weaved (fl995h and suggest the best fit 
SN yields for a-elements and iron group elements, based 
on the comparison between one zone chemical evolution 
calculation and obs ervational data. For oxy gen, they 
adopt the yield by IWooslev fc Weaver! (|1995l ) without 
modification. Their best fit yields are, as a matter 
of course, consistent with observations as far as one 
zone model with the Salpeter IMF goes. We adopt 
their yields to the hierarchical evolution model with the 
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Fig. 2. — Theoretical SNe II yields of t he zero metallic- 
ity m assive stars a gainst initial stellar mass by IKobavashi et all 
iSoO^Kreds quares mVooslev fc Weaveii 1T995D ( green filled circles) , 
IFrancois et al.l I l2004l) (blue open circles), and IChieflii fc Limongil 
l|2004l ) (magenta triangles). The open squares with dotted lines 
denote hypernovae yields. 

high mass IMF (Model KF). Since they do not discuss 
metaUicity dependence of the SN yields, the same yields 
are assumed for any metallic ity in Model KF. For Na , 
we assume the same yield as IWooslev fc Weaveil ()1995l ) 
si nce they do not de r ive be st fit one. 

IChieffi fc Limongil (j2004| ) also give explosive yields of 
massive stars from Z = to Z = Zq. Chemical evo- 
lution computation using their results is also presented 
(Model KC). We note that they chose a mass cut pa- 
rameter to eject 0.1 Mq of iron for all SNe II. Explosion 
energy is assumed to lO^lerg. 

For Type la SNe (SNe la), a single degenerate sce- 
nario is assumed. Mass range of a primaries to be SN la 
is assu med to 2Mq < m < Mup. Following the iGreggid 
(|2005l) , 9% of binaries in this mass range become SNe la. 
Delay time from star formation to SN explosion is equal 
to a lifetime of secondary companion of binary. Because 
individual stars are registered in our computation, de- 
lay time of each SN la is calculated from mass of the 
secondary star, without assuming a delay time distribu- 
tion function. Event rate and delay time distribution 
of SN la depend on the IMFs. We use the W7 model 
in llwamoto et al.l (|1999[ ) for metal yields of SNe la. In 
this study, we assume frequency and yields of SNe la are 
metallicity independent, for simplicity. It is worth not- 
ing that some studies argue that the fr action of stars to 
be SN e la depends on the metallicity. IKobavashi et al.l 
(IT998I) argue that only stars with metallicity larger than 
[Fe/H] > — 1 become SNe la. But the dependence on the 
metallicity has not yet been well understood. 

Some stars with mass around 10 Mq are thought to 
become electron capture SNe of progenitor A GB stars 
with an 0-Ne-Mg core. IWanaio et all (|2009D gives nu- 
cleosynthesys yields in this type of explosion. Amounts 
of the ejected iron and a-elements are much smaller than 
SNe II. One p r omine nt feature of the predicted yields by 
IWanaio et all (I2009D is a large yield of Zn. He argues 
that 0-Ne-Mg SNe can be a main source of Zn in the 
universe. Mass range and fraction of the stars to become 
0-Ne-Mg SNe are not well revealed. We assume that 
stars with 9 - IOMq become 0-Ne-Mg SNe. 

In this study, we omit the metal ejected by mass loss 
from intermediate massive stars. Because EMP stars arc 
formed in the very early stages of the chemical evolution, 



metal ejected from the intermediate massive stars which 
have longer lifetimes than massive stars should be negli- 
gible. Additionally, the amount of a-clemcnts, iron group 
elements, and Zn provided by the intermediate massive 
stars arc thought to be less than by SNe. 

3. OBSERVATIONAL DATA 

We adopt a bias-corrected halo MDF constructed from 
subsamples of t he HES survey with m oderate-resolution 
spectroscopy bv ISchorck et al.l jSooi). At [Fe/H] < -3, 
they have replaced the moderate resolution values with 
those derived from high resolution spectroscopy, where 
available. Since the HES survey is biased toward low 
metallicity, they evaluate biases and give the selected 
fraction, /s, as a function of metallicity and color. In- 
stead of the raw MDF, iVo,,,,([Fe/H]), of the HDS sample, 
a bias-collected MDF, iVo6s( [Fe/H])//, ([Fe/H]), is plot- 
ted to comparison with the model results. 

We also show a MDF compiled by SAGA (Stcl- 
lar Abundance for Galactic Archaeology) database 
(jSuda et al.l l2008h since we are interested in the low- 
metallicity tail of the MDF. SAGA compiles elemen- 
tal abundances of EMP halo stars which received high- 
resolution spectroscopic observations. At [Fc/H] < —3, 
metallicity derived by moderate-resolution spectroscopy 
can be significantly different with that derived by high- 
resolution spectroscopy. High-resolution data tell us ac- 
curate metallicity of EMP stars and SAGA contains a 
large sample of EMP stars. We plot raw data compiled 
by the SAGA database and this is strongly biased toward 
low metallicity. However, we may well regard the selec- 
tion of target stars for the follow-up observations as being 
hardly biased below the metallicity of [Fc/H] ~ —3. The 
MDF of Schorck ct al. (2009) contains only two stars at 
[Fe/H] < -3.6, and none at [Fe/H] < -4.3. However, 
currently 19 stars with [Fe/H] < —3.6 have been identi- 
fied in the SAGA sample and two stars have [Fe/H] < — 5. 

We consider not only the form of the MDF but also the 
total number of EMP survivors. In Paper I, we have es- 
timated efficiency of the identification of EMP survivors 
by the HES survey from its coverage area. The effec- 
tive survey area is S* = 6726 deg^ and roughly 40% of 
the candidates, selected by the objective-prism survey, 
have been examined by the spectr oscopic follow-up obser - 
vations with medium resolution (jChristlieb et al.l l20081. 
Wc assume that almost all giant EMP survivors present 
in the survey fields are detected because of the large lim- 
iting magnitude of the HES survey {B < 17.5). When 
uniform distribution of the stellar halo is assumed, ^ 5% 
of the giant EMP survivors in the Milky Way halo are 
expected to be detected. In the following figures, we plot 
the predicted number of stars which are expected to be in 
the HES sample, assuming a uniform stellar halo. When 
dc Vaucoulcurs density distribution is assumed, detection 
frequency decreases by a factor of 5, since many unde- 
tectable stars should be distributed in the inner part of 
the Galactic halo from the solar orbit. 

For elemental abundance ratio, we adopt data com- 
piled by the SAGA database (|Suda et al.ll20To| ). Figure [3] 
shows distributions of O, Na, Mg, Si, Cr, and Zn abun- 
dances relative to iron against [Fe/H]. We select high- 
resolution sample of i? > 20000. We note that, since 
SAGA compiles data from many sources, scatter of the 
abundances by the SAGA dataset can be larger than the 
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TABLE 1 
Models 



name 


IMF 


SN Yield 


KK 


Komivaet al. (2007)(M^^ = 10 Mm) 


Kobavashi et al. (2006) 


LK 


Lucatello et al. (2Q05)(M^^ = 0.79 Mm ~) 


Kobavashi et al. (2006) 


CK 


Chabrier (2003)(M^r, = 0.22 Ma^) 


Kobavashi et al. (2006) 


KW 


Komiva et al. (2007) 


Wooslev & Weaver (1995) 


KF 


Komiva et al. (2007) 


Francois et al. (2004) 


KC 


Komiva et al. (2007) 


Chieffi & Limongi (2004) 


KKn 


Komiva et al. ("2007) 


Kobavashi et al. (2006)(no hvpcrnova) 
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Fig. 3. — Abundance ratio distributions of metal poor stars compiled by the SAGA database. Among the SAGA sample, stars analysed 
by the First Stars project {blue triangles, A), stars analysed by Aoki et al., {green squares, □), and stars in Honda et al.(2004, 2007) 
{magenta inverted triangles, V) are plotted with error bars, where available. Other giants {red plus signs, +) and main sequence stars 
{yellow crosses, x) in the SAGA sample are plotted without error bars. Stars in common between two or three subsamples arc connected 
with lines. See text in detail. 

intrinsic scatter of the stellar abundances. To see the 
systematic difference between literatures, three subsam- 
ples analysed by different authors are plotted by different 
symbols with error bars. Blue triangles(A) ii i Fig . [3lde - 
note data analysed by the First Stars project (iHilll (|2002[ ) 



and other 13 papers of the series). Green squares(n) de- 
note data of which first author of source paper is W.Aoki 
(jAoki et all ()2002l ) and other 15 papers in entry list of 
SAGA). Magenta inverte d triangle (v) shows sample of 
iHonda et aiT(|2004 l2007t) . These three subsamples are 
analysed assuming a plane parallel stellar atmosphere 
(ID) and local thermal equilibrium (LTE) but there are 
systematic difference in the abundance ratio owing to 
difference in model atmospheres, parameters used dur- 
ing the analysis, and lines used in analysis. Other stars 
are plotted with red plus signs (-I-) and yellow crosses 
(x) for giants(rcff < 6000 and logg < 3.5) and dwarfs. 



respectively. 

We plot only giant stars in the following figures since 
there is systematic shift in stellar abun dances between 
giant s and dwarfs for some elements (jBonifacio et al.l 
l2009t ). Carbon enhanced stars with [C/Fe] > 0.5 are 
not plotted because their surface is thought to b e pol - 
luted by binary mass transfer. iKomiva et al.l (|2007[ ) 
show that carbon on these stars originates in the in- 
termediate massive companion stars and their surface 
abundances of O, N a and Mg are also aff ected by bi- 
nary mass transfer (jNishimura et al.l |2009[ ) . Recently, 
abmidance determinations with a non-LTE sch eme(e.g. 
[Mas hon ki na et al. 200 8; A ndrievskv et al.ll2010D or with 
3D model atmospheres(e.g.iAsplund fc Garca Predl2001l : 
iGonzlez Hernndez et all l2008| ) are carried out. Differ- 
ence of abundances determined with these models from 
ID LTE models is often of order ^ 0.5 dex. They possi- 
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Fig. 4. — Predicted MDFs assuming three different IMFs. 
The solid red, dashed green, and dotted blue lines show resultant 
MDFs of Models KK, LK, and CK, respectively. The gray his- 
togram denotesthebiascorrected observational MDF by the HES 
survev l|Schorck et al.|[2009l l. Black solid line shows the raw MDF 
compiled by the SAGA database. The predicted MDFs had been 
normalized in consideration of the survey volume o f the HES sur- 
vey. O nly the Model KK with the high mass IMF bv lKomiva et all 
11200 71 ) is consistent with the observed total number of EMP stars. 

bly diminish the shift between giant and turnoff stars and 
decrease dispersion of the abundance distribution. We 
comment about 3D / nonLTE effects for some elements in 
next section with comparison to model results. 

4. RESULTS AND DISCUSSION 
4.1. Metallicity Distribution Functions 

Figure |4] shows resultant metallicity distribution func- 
tions (MDFs) for three models using different IMFs. 
Solid red, dashed green and dotted blue lines denote 
Models KK, LK, and CK, respectively. All three models 
predict similar patterns of the MDFs but quite different 
total number of EMP survivors. This is because frac- 
tions of low-mass stars are different. As seen in Fig. |4l 
Model KK with the high mass IMF is consistent with 
observations but other models with the lower mass IMFs 
predict many more EMP survivors. We may overesti- 
mate the efficiency of the identification of EMP survivors 
by the HES survey, because we assume homogeneity of 
the Galactic stellar halo and the HES survey reaches the 
outer end of the Galactic halo. However, as seen in Fig.|4l 
the predicted number of EMP survivors for Models LK 
and CK is ~ 100 — 1000 times larger than observations 
and this large discrepancy cannot be explained by spa- 
tial inhomogeneity of the stellar halo and/or insufficiency 
of the survey. This indicates that typical mass of EMP 
population stars is significantly higher than nearby Pop. I 
stars, as shown in our earlier studies. 

Figure [5] shows dependence on the SN yields. Solid 
red, long-dashed green, short-dashed blue and dotted 
magenta lines denote results of Models KK, KW, KF, 
and KG, respectively. All model results arc similar in 
the total number of EMP survivors. 

At [Fe/H] < -3, the patterns of the MDF differ. These 
extremely metal deficient stars are very early generations 
of stars formed with metal ejected by only one or a few 
SN progcnitor(s) in their host halos. A MDF at [Fe/H] < 
—3 is sensitive to individual SN yields. The MDF of 
Model KK has a hump at [Fe/H] ~ —3.6. This is because 
iron yields of normal SNe II are tuned to 0.07 Mq in 
iKobavashi et al.l ()2006[) and metallicity of a primordial 
mini-halo with typical mass becomes [Fe/H] ~ —3.6 by 
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Fig. 5. — Same as Figure [4] but for models with different SN 
yields. The solid red, long-dashed green, short-dashed blue and 
dotted magenta lines show results of Models KK, KW, KF, and 
KG, respectively. 

a single SN. We see a smaller hump at [Fe/H] = —4 to 
—3 in other models, too. In Model KK, since energetic 
hypcrnovac blow out gas from mini-halos, the predicted 
number of stars with [Fe/H] ^ —3 to —2 is smaller than 
other models. 

Model KG shows better consiste ncy with obser v ations 
than other models. The MDF of ISchorck et al.1 (|2009l ) 
shows a steep drop around [Fc/H] = —3.6. In Model KG, 
since iron yield and explosion energy are assumed to be 
the same for all the SNe II, all the second and later gen- 
erations of stars are distributed above [Fe/H] ^ —3.6. 
However, theoretical iron yield strongly depend on the 
SN model parameters and it is difficult to distinguish 
which theoretical yield is the best one only from com- 
parison with observed the MDF. Additionally, at such 
a low metallicity, the observational sample is very small 
and the pattern of the MDF has not yet been well re- 
vealed. 

The predicted number of stars with [Fe/H] < —4.5 is 
very small and MDFs are bumpy because of numerical 
fluctuation. The number of observed HMP stars is com- 
parable to the model results. 

For more metal rich stars, a MDF depends on an aver- 
aged yield but is almost independent of characteristics of 
individual massive stars. At [Fe/H] > —2, the observed 
MDF overwhelms the predicted number of stars by the 
high mass IMF models. One possible explanation of this 
excess is a changeover of t he IMF. The IMF of m etal rich 
stars is peaks at low mass. lYamada et al.l (|2011| ) indicate 
that an IMF can be a changeover from high mass to low 
mass at [Fe/H] = —2.2 based on statistics of stellar abun- 
dance of Zn and Go. If typical mass becomes lower, the 
number of stars surviving increases. At this metallicity, 
some thick disk stars are thought to contaminate in the 
sample. It also increases the number of stars. Formation 
of the thick disk and the IMF of these more metal rich 
stars will be investigated in our future works. 

4.2. Abundance Ratio Distribution 

In Figure [6l wc show predicted and observed abun- 
dance ratio distributions for six elements for Model KK. 
The predicted typical abundance ratio and their trend 
against metallicity are si milar to the result of the chem- 
ical evolution model by IKobavashi et al.l ()2006f ). But 
since abundance ratio differs from mini-halo to mini-halo 
in our model, abundance distributions show significant 
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Fig. 6. — Distributions of elemental abundances of O, Mg, Si, Na, Cr, and Zn relative to iron against metallicity. Sampled stars in 
Model KK (red open circles) and observational sample compiled from the SAGA database (black s ymbols). Triangles (A): the First Stars 
project; squares (■): stars analysed by Aoki et al.; inverted triangles (t): IHonda et al.l I I2004I , [20071 ): plus signs (+): other giants compiled 
by the SAGA database. 



dispersion. The predicted scatters of abundances are 
prominent at [Fe/H] < —3. At larger metallicity, ele- 
ment abundances are averaged by mixing of ejecta from 
many SNe and the scatters decrease. As shown in Pa- 
per I. metals of stars with [Fe/H] = —2.5 originate in 
^ 10 precursory SNe in their host mini-halos. 

Most stars with [Fe/H] < —4 are local first stars with- 
out any precursory SN exploded in their host mini-halos. 
Metals in these stars are originate in mixtures of mat- 
ter ejected to IGM by SNe which had exploded in other 
mini-halos prior to the formation of the host mini-halos 
of these stars. Scatter of the abundances of these stars 
is smaller than the stars with —3 > [Fe/H] > —4. Ob- 
served HMP stars show abundance anomalies of C, N, 
O and Na but these are thought to b e due to the sur- 
face pollution by bina ry mass transfer (|Suda et al.ll2004l : 
iNishimura et al.ll2009| ). 

4.2.1. a-elements 

For Mg and Si, typical abundance ratio is consistent 
with observations. Observed distribution of [a/Fc] is al- 
most flat against [Fe/H] at [Fe/H] < — 1 and the disper- 
sion is small. The predicted dispersion of the a-element 
abundances is comparable to or smaller than the disper- 
sion of the observational sample. 

Observationally, [0/Fc] of EMP stars shows a slight 
increasing trend as metallicity decreases but predicted 
abundances do not. Typical O abundance of the First 
Stars sample is ^ 0.2 dex higher than the predicted 
value. As seen in Fig. [^ a larger amount of O is yielded 
in a more massive progenitor. If this increasing trend 



is real, stars more massive than 20 Af© thought to be 
the dominant source of metal in EMP stars. However, 
Mg and Si abundances show no such a trend although a 
larger amount of these elements is yielded in a more mas- 
sive progenitor, too. Another possible additional source 
of oxygen for EMP stars is intermediate massive AGB 
stars. But AGB stars eject both O and C, O abundances 
of stars other than CEMP stars should not strongly af- 
fected by AGB stars. We note that observational deter- 
mination of the oxygen abundance is diffic ult and there 
is larg e uncertainty on observational data. iNissen et al.l 
(l200l argue that, when 3D effects are taken into ac- 
count, [0/Fe] decreases in the metal poor stars and the 
increasing trend diminishes. For some EMP stars with 
low [0/Fc], oxygen is undetectable. 

The sample from the First Stars project shows very 
small scatter for Mg and Si. This is consistent with the 
result of Model KK. As seen in the following. Model KK 
with hypernovae predicts the smallest scatters among 
the computations in this paper. It is because [0/Fe] 
in the hypernovae ejecta is tuned to be constant and 
large energy of the hypernovae efficiently mix the ISM. 
The observed small scatter of the First Stars sample in- 
dicates that gas in the early universe is mixed with large 
mass. However, the subsamples by other authors show 
larger scatter (~ 0.5 dex). Further observations are re- 
quired to understand the gas dynamics in the early uni- 
verse. For Mg, the First Stars sample are distributed at a 
slightly lower [Mg/Fe] t han the model results bu t other 
samples are consistent. lAndrievskv et al.l (|2010l ) argue 
that, when non-LTE effects are taken into account, the 
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mean value of [Mg/Fe] increases and become similar to 
the mean value of [0/Fe]. This will improve the consis- 
tency between model and observed data. 

Observationally, it is known that [a/Fe] decreases as 
[Fe/H] increases at [Fe/H] > —1 for stars in the Milky 
Way. In the many chemical evolution studies, this de- 
creasing trend_hBS been explained by the contribution of 
SNe la (jMatteuccill2008L and references therein). SNc la 
eject large amounts of iron after long delay time and de- 
crease [a/Fe] at large metallicity. Under the high mass 
IMF, the contribution from SNe la is relatively weak be- 
cause the number of intermediate massive stars relative 
to high mass stars is small. In Model KK, [a/Fe] de- 
creases at [Fe/H] > —1 but the rate of the decline is 
smaller than observations. Although the IMF is assumed 
to be the same for the whole metallicity range for sim- 
plicity in this study, our high mass IMF is derived from 
the statistics of the EMP survivors and the IMF should 
be different at higher metallicity. 

4.2.2. Na 

For Na, observational sample of EMP stars show quite 
large scatter. All the subsamples shows large scatter 
and it is thought to be intrinsic. Observed stars with 
lower [Na/Fe] have abundance similar to the model re- 
sult but some other stars show larger [Na/Fe]. This in- 
dicates that there are additional sources of Na. Internal 
mixing may modify the surface Na abun dance of some 
evolved EMP survivors jSpite et al.ll2006D . As discussed 
in lNishimura et ahl (|2009[ ). in intermediate massive EMP 
stars at AGB stage, Na can be synthesized and dredged- 
up by the He-flash driven deep-mixing, which is a mixing 
mechanism peculiar to the EMP stars. Stars with high 
[Na/Fe] can be influenced by matter ejected from these 
AGB stars. At [Fe/H] > —2, predicted typical abun- 
dance and their increasing trend are consistent with ob- 
servations. We note that, when non-LTE effect taken 
into account, Na abundance of some stars for which high 
[Na/Fe] value is reporte d decreases to [Na/Fe] ^ —0.2 
i|Andrievskv et al.l[2007l ). 

4.2.3. Cr 

[Cr/Fe] shows a decreasing trend as decreasing metal- 
licity. The predicted relative abundance is not consis- 
te nt with the observations of EMP stars, as pointed out 
bv iKobavashi et al.l ()2006l ). We note that the sample of 
iHonda et all (|2004[ ) is distributed around [Cr/Fe] = 
and consistent with model results. However, the First 
Stars sample and Aoki's sample show clear decreasing 
trend with small scatter. 

4.2.4. Zn 

In the models in this paper, Zn is r nainly produced 
by 0-Ne-Mg SNe since IWanaio et all (|2009l ) predict a 
large Zn yield for 0-Ne-Mg SNe. Hypernovae are also 
important sources of Zn, especially at very low metallic- 
ity. For EMP stars, all observational subsamples show 
similar abundance distribution and they are consistent 
with the model result. 

Observations show an increasing trend of [Zn/Fe] as de- 
creasing metallicity. Both subsamples by the First Stars 
group and Aoki et al. also show a clear increasing trend 
with small scatter. The model result shows flat distri- 



bution and predicts lower [Z n/Fe] at [Fe/H] > —2. Re- 
centry, lYamada et al.l (j2011| ) shows that a decrease in 
[Zn/Fe] above [Fe/H] > —2.2 may be due to changeover 
of the IMF from high mass to low mass. The IMF 
changeover lowers the frequency of hypernova and low- 
ers the [Zn/Fe]. The contribution from hypernovae is 
discussed again with results of the model without hyper- 
novae in Section l4?3l 

At [Fe/H] = —2 to —1, the model predict higher 
[Zn/Fe] than observations. Zn thought to be overpro- 
duced by 0-Be-Mg SNe. Criteri on to be 0-Ne-Mg SNe 
is not yet revealed (jHerwig||2005| ) and the number of O- 
Ne-Mg SNe can be smaller. 

Predicted scatter of Zn is larger than other elements. 
This is because hypernovae and 0-Ne-Mg SNe eject 
matter with high [Zn/Fe] but normal SNe with mass 
m > 20 Mq eject a small amount of Zn. 

4.2.5. IMF dependence 

Figures [7] and [5] show predicted abundance ratio dis- 
tributions by models using different IMFs. Model CK 
predicts similar typical abundance ratio to Model KK. 
O, Mg and Si, are mainly provided by stars with mass 
heavier than 20 Mq and iron is ejected by all SNe II and 
SNe la. For EMP stars, [a/Fe] depends on a fraction 
of stars with > 20 Mq among SNe II and it depends 
on the IMF. For Models KK and CK, typical mass of 
stars is quite different but the slope of the IMFs at mass 
range to be SNe II (10 — 50 Mq) is similar, as seen in 
Fig. [TJ Relative frequency of heavier (> 20 Mq) and 
hghter (< 20 Mq) SNe II are similar for both IMFs 
and typical [a/Fe] are also similar. However, Model 
CK predicts that some stars with low oxygen abundance 
([0/Fe] ~ 0.2 — 0.3) are distributed at very low metal- 
licity range ([Fe/H] < —3.5). These stars are born in 
mini halos formed at low redshift. In this model, many 
SNe la yield iron at lower redshift, and iron ejected from 
mini halos lower the [0/Fe] of IGM. Metalhcity of the 
IGM is still low because ejected matter is diluted in 
large mass. Mini halos formed with the IGM polluted 
by SNe la have low [O /Fe] but low metallicity. Observa- 
tionally, these stars are not detected. If ejected matter 
is mixed in smaller mass, iron and oxygen abundance of 
polluted IGM become larger and these 0-poor stars dis- 
sipate. We note that, oxygen abundance of these stars 
can be lower than the detection limit of O. Some stars 
without detection of O possibly have such an abundance 
feature. 

For Model LK, [0/Fe] is obviously lower than Models 
KK and CK in the whole metallicity range, as seen in 
Fig. |5| This is because a slope of the IMF is steeper at 
mass range of stars to be SNe II. The relative frequency 
of heavier SNc II (> 20 M©) is smaller and a smaller 
amount of a-elements is ejected. The predicted distri- 
bution of [O /Fe] is much lower than the observations for 
EMP stars. At [Fe/H] ~ -2, the abundances of Mg, Si 
and Na relative to iron are also lower than the observa- 
tions. 

At [Fe/H] > —1, we can see clear decreasing trends for 
the a-element abundances as increasing metallicity for 
Model CK and LK. Relative numbers of SNe la are larger 
for these models and they lower the [a/Fe] at higher 
metallicity. 
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4.2.6. Supernova Models 

Figures m Uni and [TT] show abundance ratio distri- 
bu tions for Models KW, KF and KC using SN y ields 
bvlWooslev fc Weaveij(ll995D .lFrancois et al.l (l2004f). and 



IChiefR fc Limongir (2( i04), respective ly. 

For Model KW, as [Francois et al.l (I2004D pointed out, 
predicted a-element abundances do not agree with ob- 
servations. Predicted abundances of the O, Mg, and Na 
relative to iron are 0.5dex or more lower than the obser- 
vational sample. 

For Model KF, since [Francois et ahl ()2004D modify the 
nucleosynthetic yields to match observational data, pre- 
dicted typical abundances of EMP stars show good agree- 
ment with observations for elements other than Na. The 
decreasing trend of [Cr/Fe] at low metallicity is also re- 
produced. They assume that a SN with larger initial 
mass ejects a smaller amount of Cr (see Fig . Since a 
star with larger mass have shorter lifetime, [Cr/Fe] in- 
creases as metallicity increases with time. However, all 
other studies with nucleosynthesis computations predict 
that a more massive star yields a larger amount of Cr. 

While Model KF well reproduces the typical abun- 
dances, it predicts larger dispersion of the element abun- 
dances than the observations. Especially, [0/Fe] dis- 
tributes from -1.5 to -1-1.5 at [Fe/H] < -3. This model 
also predicts some stars with [Si/Fe] > +1. The pre- 
dicted scatter is much larger than the observations and 
it indicates that a one-zone model is inadequate to un- 
derstand earliest phases of the chemical evolution and 
metal yields of very metal poor SNe. 

Model KC predicts typical abundances in agreement 
with observational sample, for Mg and Si. Lower O 
abundance than observations is predicted but observed 
O abundances can be lower when non-LTE and 3D effect 
taken into account, as mentioned above. 

This model predict some stars with very low [0/Fe] 
and [Mg/Fe] at [Fe/H] < —3. Such abundance pat- 
terns are produced from SNe at low-mass end of the 
mass range to be SNe II. As seen in Figure 2, stars with 
10 — 12 Mq yield small amount of O and Mg. Lower mass 
limit to be SNe II is assumed to be 10 Mq in this paper 
but fate of the stars with ~ 10 Mq is not well revealed. 
Some stars wi th '-^8—12 Mp, thought to become "super- 
AGB" stars (jGarcia-Berro fc IbeiJll994f ) and evolve to 
O-Ne-Mg whi te dwarfs or electron capture supernovae 
(jHerwigl l2005f ) with very little iron yield. Absence of 
the very a-poor stars possibly indicate that the lower 
mass limit to be SNe II is larger tha n 10 Mq at very 
low metallicity. iKawabata et al.l ()2005() argue that stars 
with 8 — 12 Mq become "faint supe r novae " with low iron 
yield. Although IChiefE fc Limongil (|2004[ ) have assumed 
iron yield is 0.1 M0 for all SNe, observations indicate 
that some SNe yield lower amounts of iron. 

Na is overproduced at higher metallicity ([Fe/H] > 
—2). Cr abundance at solar metallicity is consistent with 
the observations but the increasing trend is not repro- 
duced. 

Very large scatter of Zn is predicted because they argue 
that a large amount of Zn is yielded in a SN II with 
m < 13 Mq but a very little amount of Zn is yielded in a 
SN with TO > 13 Mq. When low mass limit to be SNe II 
is larger as discussed above, stars with very high [Zn/Fe] 
is not formed. A yield of Zn is sensitive to entropy during 



explosive Si-burning at SN explosion. As discussed later, 
energetic hypernovac thought to be required to explain 
Zn abundance of EMP stars and their trend. 

4.3. Hypernova v.s. normal SN 

Figures [12] and [13] show a result of Model KKn with- 
out hypernovac contribution. The MDF of Model KKn 
is similar to Model KK. For the abundance ratio distri- 
butions, the most plausible difference from Model KK 
is lower [Zn/Fe]. Since hypernovae synthesize a much 
larger amount of Zn than normal SNe, Zn abundance 
of Model KKn is lower than Model KK and lower than 
observations for EMP stars, a-clcmcnt abundances pre- 
dicted by Model KKn are slightly higher than Model 
KK and [0/Fe] shows better agreement with observa- 
tions than Model KK. But for Mg and Cr, Model KKn 
predicts slightly higher relative abundances than the ob- 
servations. This is because a normal SN yields a smaller 
amount of iron than a hypernova. 

Large explosion energy of hypernovae affects also gas 
dynamics. Many mini-halos are blown up by their large 
explosion energy and ejected metal is mixed in a large 
mass. Since it averages element abundances, the scat- 
ter of the predicted abundance of Model KK is smaller 
than other models. The observed small scatter of [a/Fe] 
suggest that there were many hypernovae in the early 
phas es of the chemica l evolu tion. We note that, how- 
ever, [Kobavashi et al.[ (]2P06!) have tuned parameters in 
their computations to get [0/Fe] = 0.5 for all hypernovae 
and the scatter of the a-elemcnt abundances is decreased 
artificially. 

5. CONCLUSIONS 

We compute formation history of EMP stars with a hi- 
erarchical chemical evolution model and present MDFs 
and abundance ratio distributions. We adopt vari- 
ous IMFs and SN yields and compare the results. As 
for IMFs, we adopt the hig h mass IMF given in our 
earher studies ([Komiva et al.l [20071 [2009aD . the IMF by 
[Lucatello_et al.[ ([2005[) . and the standard low mass IMF 
by [Chabriej ([20031). The pattern of the MDFs is simi- 
lar for these three models but the predicted numbers of 
EMP survivors are quite different. The high mass IMF 
model predicts a comparable number of EMP survivors 
to observations but other two IMFs predict many more 
EMP stars in our scenario. Our hierarchical model re- 
produce a steep dechne of the MDF below [Fe/H] = —3.6 
and a tail below [Fe/H] < -4. 

Abundance ratio distributions predicted by the high 
mass IMF are similar to the Salpeter IMF and the pre- 
dicted a-element abundances are consistent with obser- 
vations. Typical value of [a/Fe] depends on the slope 
of the IMF at a mass range in which stars explode as 
SNe II and the high mass IMF with Mmd = IOA/q and 
Am = 0.4 has a slope s i milar to the Salpeter IMF. The 
IMF bv [Lucatello et al] (|2005[ ) which has a steeper slope 
at 10 — 50 M0 predicts lower [a/Fe]. 

We investigate not only the averaged abundance and 
their trend against metallicity but also the scatter of the 
abundances using our hierarchical models. Abundance 
distributions strongly depend on the adopted SNe nucle- 
osynthesi s models. For typ i cal ab u ndances of a-el e ment s , 
yields by 'Kobavashi et al.' ('2 006D . [Francois et al.l ()2004D . 
and [Chicffi fc Limonai (200^ show reasonable agree- 
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Fig. 12. — Effect of hypcrnovae on the MDF. The solid red 
and dashed blue lines denote result of Model KK and Model KKn 
without hypernovae, respectively. 

ment with the observational sample. The sample from 
the First Stars project shows very small scatter for 
[Mg/Fc] and [Si/Fe]. This indicates that abundance ra- 
tio ejected by various SNe is homogeneous or inter stel- 
lar matte r is well mixed in large mass. The result with 
yields by iKobavashi et all ()2006[ ) is consistent with this 
very small scatter. However, stars analysed by other au- 
thors show larger sc atter. Although the SN yields by 
[Francois et al.l ()2004l ) are the best fit yields as far as aver- 
aged abundances, their yields predict much larger scatter 
than all the observational sample for O and Si. Absence 
of stars with low [0/Fc] and low [Mg/Fe] possibly indi- 
cate that iron yields of stars with mass around 10 M© 
are lower than normal SNe H. 



Observed decreasing trend of [Cr/Fc] as mctallicity de- 
creases cannot be explained by the adopted yields by the 
theoretical SN nucleosynthesis computations. For Na, 
intermediate massive AGB stars or internal mixing af- 
fect surface abundances of some EMP stars with Na en- 
hancement. For O, Mg and Na, correction of 3D and/or 
non-LTE effect thought to decrease the scatter of ob- 
served abundances and improve the consistency to the 
model result. Hypernovae are the plausible dominant 
source of Zn for EMP stars. Models without hypernovae 
predict lower [Zn/Fe] than in observations. Observed in- 
creasing trend of [Zn/Fe] as metallicity decreases is not 
reproduced in our hierarchical models. 

Results in this paper are depend on the models of 
SN nucleosynthesis and accuracy of the abundance de- 
termination from spectroscopic observational data. Un- 
certainty of the theoretical SN yields is still large. Ob- 
servational subsamples by different authors show differ- 
ent abundance distributions for some elements. Recent 
studies with 3D/non-LTE effects taken into account show 
that averaged abundances shift significantly by these ef- 
fects and dispersion of abundances become smaller for 
some elements. Further theoretical and observational 
studies are required to reveal the origin of metals in the 
EMP stars and nature of the metal poor SNe. 



APPENDIX 
PAIR INSTABILITY SUPERNOVA 

First stars formed without metal arc thought to be very massive stars with mass larger than 100 Mq. Many numerical 
simu lations for first star for ming cloud s argue that onl y one very mas sive density peak is formed without fragmentation 
(e.g. lAbel et al.. 2002: O'Shea fc Norman 2007: Yoshida et al.ll2008l but see also Clark et al. 2008). Stars with mass 
140 — 270 Mq become PISNe with huge explosion energy. Theoretical studies of SN nucleosynthesi s show that metal 
yields by PISNe are quite different from SNe II (|Umeda fc Nomotol[200l iHeger fc Wooslevll2002[ ). The amount of 
metal ejected by the energetic PISNe is much larger than from SNe II. They show a strong "odd-even effect"; i.e., 
ratio between the yields of elements with even atomic number and odd atomic number is quite large. Zn production 
is very inefficient relative to other elements such as iron. However, observational studies of the element abundances of 
the EMP stars show that there is no nucleosynthetic signature of the PISNe. 

Metal provided by massive stars is thought to change the IMF and enable low mass star formation. The critical 
metallicity to the low mass star formation is investigated through studies of cooling processes in molecular cloud by 
dust and/or metal. Studies with dust argue that cooling by H2 molecule fo rmed on the surface of the du st sufficiently 
cool the gas cloud and lower their Jeans mass at metallicity [Fe/H] > —6 (jOmukai et al.|[2"005l : [Schneid er 2006), and 
low mass stars are formed. 

In this appendix, we present a computation with PISNe taken into account. We assume that typical mass is very 
massive, Mmd = 200 M©, for stars with [Z/H] < -6 an d stars with mass 140 - 270M(d become PISNe. Other 
assumptions are the same as for Model KK. PISN yields bv lUmeda k Nomotol (|2002[ ) are adopted. 

Figures ET) and K2] show MDFs and relative abundance distributions, respectively. The results arc quite similar to 
the results of Model KK without PISNe and we cannot distinguish the observable signature of the PISNe. PISNe eject 
a large amount of metal at energetic explosions and the ejected metal is mixed into IGM. Metallicity of the IGM is 
enriched over [Z/H] = —6 by a few dozen PISNe occurring at very high redshift. After that, the IMF is changed to the 
high-mass one for EMP stars (Mmd = 10 Mq). Matter ejected by SNe II is mixed into the IGM and SNe II become the 
dominant source of the metal for stars with [Fe/H] > —5. Additionally, as we have shown in Paper I, surface pollution 
by accretion of interstellar matter changes surface abundances of HMP stars with [Fc/H] < —5 after their formation 
and the signature of the PISNe should be covered. In conclusion, the abundance pattern peculiar to the PISNe ejecta 
is obscured. 
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